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1.  INTRODUCTION 


The  work  reported  here  was  done  at  the  Ballistic  Modeling  Division  of  the  Ballistic  Research 
Laboratory,  (BRL)  in  support  of  theoretical  and  experimental  efforts  througout  the  BRL.  The 
objective  of  this  work  was  the  refinement  and  maintenance  of  the  Eulerian,  hydrodynamics 
code,  BRLHELP1,  at  the  state  of  the  an,  for  application  to  problems  in  both  basic  research  and 
experimental  studies.  This  objective  was  approached  through  the  analysis  of  calculations  of 
problems  of  current  interest,  and  of  the  applicability  of  various  models  in  BRLHELP  to  these 
problems. 

The  HELP2  code,  an  evolutionary  antecedent  of  BRLHELP,  has  been  used,  for  some  time, 
as  a  predictive  tool  in  warhead  mechanics,  and  intermediate  and  terminal  ballistics.  Tts  use  is 
limited,  however,  by  serious  problems  with  the  accuracy  and  robustness  of  some  of  its  models. 
BRLHELP  is  devoid  of  most,  if  not  all,  of  these  problems.  It  has  been  modified  to  calculate 
internal  energy  accurately3,  has  more  robust  treatments  of  certain  numerical  problems  inherent 
in  the  old  HELP  code,  and  has  temperature  available  as  an  output  variable  for  metals  and  alloys 
via  the  BRLGRAY4  equation  of  state. 


2.  BRIEF  DESCRIPTION  OF  MODELS 

Detailed  descriptions  of  the  mathematical  and  numerical  models  in  BRLHELP  can  be  found 
in  References  1-3.  A  brief  description  is  given  here. 

BRLHELP  is  a  two-dimensional,  Eulerian  hydrodynamics  code.  It  describes  unsteady 
multi-material  problems,  treats  material  strength  as  an  elastic-plastic  phenomenon,  and 
incorporates  a  multiphase  EOS  for  metals. 

The  code  uses  massless  tracer  panicles  to  define  material  interfaces  and  free  surfaces.  These 
particles  are  placed  along  surface  boundaries  and  move  with  the  local  material  velocity.  This 
Lagrangian-like  description  provides  an  accurate  definition  of  the  material  interfaces  and  free 
surfaces  without  sacrificing  the  ease  of  treatment  of  extreme  distortion  by  the  Eulerian  model. 

Input  to  the  code  consists  of  an  initial  geometric  configuration,  initial  values,  explosive 
loading  parameters,  material  properties,  and  desired  editing  features. 

The  geometry  package  can  describe  any  configuration  decomposable  into  straight  lines,  arcs 
of  circles,  and  arcs  of  ellipses.  The  high-explosive  detonation  model  can  treat  multiple 
explosives,  multiple  detonation  points,  and  the  presence  of  wave  shaping  devices.  Sliplines  may 
be  employed  at  material  interfaces  to  allow  different  materials  to  develop  velocities 
commensurate  with  their  material  properties,  rather  than  a  common  interface  velocity. 


1.  Joseph  Lacetera ,  Janet  E.  Lacetera,  James  A.  Schmitt ,  'The  BRL  7600  Version  of  the  HELP  Cod?  Ballistic  Research  La¬ 
boratories  Report:  ARBRL-TR-02209,  Jan  SO.  (AD  #A082804) 

2.  Laura  J.  Ha&man,  et  al "HELP,  A  Multi-material  Eulerian  Program  for  Compressible  Fluid  and  Elastic-Plastic  Flows  in 
Two  Space  Dimensions  and  Timf,  Systems ,  Science ,  and  Software  Report  TR-76-45-BK2 ,  Apr  1976. 

3.  James  A L  Schmitt ,  Truncation  Error  Terms  in  the  Kinetic  Energy  Calculation  in  the  HELP  Algorithm  and  Their  Conse¬ 
quence^,  J.  Comp.  Phys,  Vol  35,  No  2.  Apr  1980 . 

4.  Joseph  Lacetera  "BRLGRAY:  The  Ballistic  Research  Laboratory  Version  of  the  GRAY  Equation  qfStatf,  ARBRL-TR- 
02258,  Aug  1980,  (AD  #A091246) 
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Output  includes  pressure,  density,  material  velocity  vectors,  specific  internal  energy,  and 
deviatoric  stress  components  as  functions  of  space  and  time.  Temperature  of  metals  and  alloys 
is  also  available  when  BRLGRAY  is  chosen  as  the  equation  of  state. 


3.  APPLICATIONS 

BRLHELP  is  a  necessary  numerical  tool  for  the  performance  of  optimization  and 
prefabrication  studies  for  materials  and  geometries  proposed  in  shaped-charge  warhead  design. 
It  has  the  capability  of  modelling  the  collapse  of  the  liner  and  the  formation  and  development 
of  the  jet  and  the  slug.  It  provides  information  which  can  be  used  to  determine  jet-tip  velocity, 
collapse  velocity  and  collapse  angle,  liner  position  as  a  function  of  time,  stretching  rate  of  the 
jet,  and  strength  effects. 

The  HELP  code  had  been  used  with  success  in  predicting  jet-tip  velocities  and  liner  collapse 
for  a  variety  of  geometric  designs5-7.  However,  HELP  overestimated  the  internal  energy  in  the 
collapsed  portion  of  conical  shaped-charge  liner  materials.  This  overestimation  led,  in  some 
cases,  to  unreasonably  low  estimates  of  jet  densities5.  Because  of  this  an  artificial  bound  on  the 
internal  energy  in  metals  in  shaped  charge  problems  was  installed  in  the  Tillotson8  EOS  in  an 
attempt  to  maintain  a  more  appropriate  jet  density.  Without  this  bound,  the  jet  tip  exhibited 
internal  energies  so  high  and  densities  so  low  as  to  imply  vaporization. 

In  BRLHELP,  these  problems  have  been  alleviated.  Internal  energies  are  calculated  quite 
accurately,  and  the  artificial  bound  on  internal  energy  in  metals  in  shaped-charge  liner  problems 
has  been  removed.  Calculations  are  performed  without  artificial  constructs  and  physically 
reasonable  results  are  being  obtained,  including  experimentally  verifiable  metal  temperatures. 

4.  UPDATED  VERSIONS  OF  THE  BRLHELP  CODE 

All  of  the  input  to  BRLHELP,  except  for  an  alphanumeric  header  card,  has  been  converted 
to  the  NAMELIST  format. 

The  NAMELIST  input  streams  are  called: 

OPTNS,  START,  RUN,  GRID,  MTRL,  GEQM,  SLIP,  and  DETN. 

In  the  HELP  code  input  variables  located  in  the  first  150  words  of  blank  common  are  called 
Z-variables  because  they  follow  Z(l),  the  first  word  of  blank  common,  and  therefore  can  be 
referenced  as  Z(T).  A  set  of  control  parameters  are  first  read  using  the  NAMELIST:  OPTNS. 
The  Z-variables  are  then  read  using  the  NAMELTSTs:  START  and  RUN. 


5.  J.  T.  Harrison,  " Comparison  Between  the  Eiderian  Hydrodynamic  Computer  Code  (BRLSC)  and  Experimental  Collapse  of  a 
Shaped  Charge  Li  net'',  BRL  Memorandum  Report  ARBRL-MR-02841 ,  Jun  1978  (AD  #059711). 

6.  Janet  E.  Lacetera  and  William  P.  Walters,  'Theoretical  and  Experimental  Studies  of  Hemispherical  Shaped-Charge  Liner. f. 
Transactions  of  the  25th  Conference  of  Army  Mathematicians,  ARO  Report  80-1,  Jan  1980. 

7.  A.  Kiwan,  and  A.  Arbuckle.  "Study  of  Liner  Collapse,  Jet  Formation  and  Characteristics  from  Impulsive  ShapedCharge  Sys- 
■  temf.  BRL  Report  No.  2028,  Nov  1977  (AD  #051342). 

8.  J.  H.  Tillotson,  "Metallic  Equation  of  State  for  Hypervelocity  Impact f,  General  Atomic  Report  No.  GS-3216,  Jul  1962. 
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The  fatal  error  caused  by  pure  cells  overemptying  has  been  corrected  using  coding 
generated  by  Janet  Lacetera9.  This  coding  is  discussed  in  Section  6  of  this  report  Pressure 
iteration  problems,  inherent  in  old  versions  of  the  code,  are  discussed  in  Section  7  and  a  new 
mass-weighted  average  pressure  scheme  is  discussed  in  Section  8.  A  call  to  a  subroutine  called 
BRLEOS  has  been  installed  in  the  old  subroutine  EQST,  in  place  of  equation-of-state  coding,  to 
allow  access  to  the  HRLEOS  code  system10  which  includes  BRLGRAY,  the  three-phase  EOS  for 
metals. 


There  have  been  several  versions  of  BRLHELP  active  on  the  BRL  7600.  They  are  cycles 
14  through  23  of  the  permanent  file  BRLHELP,  with  ID=CMCHYDRO.  They  differ  as 
follows: 


CY-14 
CY-15 
CY=16 
CY=17 
CY=18 
CY=19 
CY=20 
CY-21 
CY=22 
CY  “23 


BRL  7600  version  of  UNIV AC-1 108  HELP2. 

Cycle  14  modified  for  accurate  energy  calculation3. 

Cycle  14  with  complete  NAMELIST  input 
Cycle  15  with  complete  NAMELIST  input 

Cycle  17  with  OVERMT  correction  and  temporaxy  pressure  iteration  correction. 

Cycle  18  with  FTN  compiler  modification 

Cycle  19  with  mass- weighted  average  pressure  for  mixed  cells. 

Cycle  20  with  most  of  EQST  functions  transferred  to  BRLEOS 
Cycle  21  with  remaining  EQST  functions  added  to  BRLEOS. 

Cycle  22  with  temperature  map  for  metals. 


Cycles  14  through  22  have  been  retired. 


A  listing  of  the  jobstream  which  created  cycle  17  of  BRLHELP,  after  cycle  16  had  been 
created,  is  shown  in  Reference  1.  Listings  showing  the  jobstreams  for  creating  later  cycles  are 
shown  in  Appendix  A  of  this  report. 


5.  DISCUSSION  OF  INPUT 

A.  Capabilities  and  Limitations 

The  problem  generator  allows  the  user  to  generate  material  interfaces  which  are  composed 
of  straight  line  segments,  or  portions  of  circles  or  ellipses,  with  the  following  limitations: 

1.  The  number  of  material  packages  in  a  calculation  is  limited  only  by  the  dimensions  of  the 
material  arrays. 

2.  A  single  cell  can  contain  any  number  of  material  interfaces  and  any  number  of  distinct 
materials. 

3.  Different  material  packages  can  contain  the  same  material  with  different  initial  conditions 
(velocities,  energies,  densities). 

4.  A  package  can  be  divided  into  any  number  of  spatially  distinct  subpackages. 


9.  Janet  Lacetera,  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground  Private  communication 

10.  Joseph  Lacetera  and  William  Buchheister,  * BRLPUFF  and  the  BRLEOS  Code  Systenf,  Ballistic  Research  Laboratory  In¬ 
terim  Memorandum  Report  No.  655,  Aug  1979. 
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5.  The  material  in  each  package  must  be  homogeneous. 

6.  The  grid  boundaries  do  not  act  as  free  surfaces.  The  code  assumes  that  a  package  that 
extends  to  the  edge  of  the  grid  is  infinite  in  that  direction 

7.  The  grid  must  have  at  least  three  rows  and/or  columns. 

8.  A  calculation  can  be  made  in  cylindrical  or  rectangular  coordinates. 

9.  The  top  and  right  grid  boundaries  are  transmittive.  The  left  boundary  is  reflective  (z-axis 
in  cylindrical  coordinates).  The  bottom  boundary  can  be  either  transmittive  or  reflective. 


B.  Edit  Control  Parameters  and  Z-Variables 

Edit  control  parameters  are  read  using  the  NAMELTST:  OPTNS.  The  Z-variables  are  then 
read  using  the  NAMELISTS:  START  and  RUN.  These  are  defined  below. 

NAMELIST  OPTNS:  (Edit  control  parameters;  it  should  be  noted  that 
the  first  three  of  these  parameters  suppress  diagnostic  print  and 
should  be  used  with  caution) 


TFLGST: 


TPMADJ: 


—0,  default 

— 1,  suppresses  message  "SUM  OF  FRACS  NOT  EQUAL  TO  TOTAL..." 
— 0,  default 

— 1,  suppresses  message  "EVACUATION  OF  MATERIAL  N..." 


ITPHSR 


— 0,  default 

—  1,  suppresses  message  "MASS  EVAPORATED  DUE  TO  ROUNDOFF..." 


NPRTOP: 


“0,  default 

“1,  suppresses  edit  of  sliplines 
-=2,  suppresses  long  edit  on  last  cycle 
—3,  combines  effects  of  (1 )  and  (2). 


TTMMAX:  The  maximum  time  in  CPU  seconds  (decimal)  that  the  code  will 

use  in  the  calculation 

For  NPRTOP  —2  or  3  TTMMAX  should  be  about  4  seconds  less 
than  the  jobcard  time;  otherwise  TIMMAX  should  be  about  10 
seconds  less  than  the  jobcard  time. 
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NAMELIST  START 


PK(1): 

Problem  identifier  (numeric) 

PK(2): 

Restart  cycle  number 

PK(3): 

-0 

for  startup  run; 

—1 

for  restart  without  short  edit  print  of  restart  cycle; 

—2 

for  restart  with  short  edit  print  of  restart  cycle. 

NAMELIST  RUN 

Contains  all  Z- variables  not  contained  in  NAMELIST  START. 


C.  Cell  Dimensions 

Input  data  which  define  cell  dimensions  are  read  into  NAMELIST  GRID.  Dimensions  of 
cells,  in  the  area  of  interest,  should  not  vary  by  more  than  10%,  nor  should  large  aspect  ratios 
(e.g.  dx/dy  >  2)  be  used  near  the  region  of  interest 

D. 

Initial  Conditions 

In  order  for  the  material  tracer  particle  routine  to  function,  the  flag  of  the  cell  in  the  bottom 
left  corner  of  the  grid,  MFLAGK2),  and  the  material  code  number  must  be  defined.  If  the  cell 
is  a  pure  cell  of  material  N,  set  MFLAG(2)—N;  if  it  will  be  an  interface,  then  set 
MFLAG(2)=-1;  if  it  will  be  a  void  cell,  set  MFLAG(2)!=0.  (In  the  case  of  MFLAG(2)=-1  its 
value  will  be  set  >  100  after  the  interfaces  are  processed  by  CAUFRC  and  VOLFND.) 

MFLACX2),  material  code  numbers  (listed  in  EQST),  internal  energy  (ergs/gm),  and  radial 
and  axial  velocities  (cm/s)  are  all  read  into  NAMELIST  MTRL.  The  equation- of -state 
constants  are  defined  in  EQST  and  COMDIM.  Variations  to  the  initial  conditions  would  require 
modifications  to  FTLGRD.  The  initial  density  RHOIN  may  be  different  than  the  material’s 
normal  density  RHOZ,  defined  in  COMDIM. 


E.  Strength  Constants 

Shear  yield  strength  and  tensile  strength  of  each  material  package  are  also  read  into 
NAMELTST  MTRL.  If  a  material  has  no  strength,  AMDM(N)  must  be  set  to  1;  the  shear 
yield  strength  constants  will  be  set  to  zero  automatically  if  not  read  into  MTRL. 
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F.  Material  Tracer  Particles 


Material  package  boundaries  are  described  by  a  series  of  straight  lines  and/or  arcs  of  cirles, 
and/or  arcs  of  ellipses.  Input  is  read  by  TSETUP  into  NAMELIST  GEOM,  and  is  defined 
below.  Data  is  read  into  GEOM  for  each  segment,  and  is  read  in  an  order  that  places  the 
package  to  the  left  of  the  segment.  Tracer  particles  are  not  placed  along  package  boundaries 
that  are  also  grid  boundaries. 


LT: 


—  1  straight  horizontal  line; 

—  2  straighi  vertical  line; 

—  3  straight  diagonal  line; 

—  4  arc  of  circle 

—  5  arc  of  ellipse 


(LT  must  be  negative  for  last  segment  of  a  package  or  subpackage; 
LT  must  equal  100  after  all  boundaries  have  been  defined.) 


MPN:  Material  package  number  (void  —  NMAT  +  1) 

NPTS:  Number  of  points  to  be  placed  along  line  or  arc 


XI S:  starting  X-coordinate 

Y1S:  starting  Y-coordinate 

X2F:  final  X-coordinate 

Y2f:  final  Y-coordinate 

TH1:  starting  angle  of  arc 

TH2:  terminal  angle  of  arc 

SLA:  X-coordinate  of  center  of  circle;  or  eccentricity  parameter  A 

in  equation  of  ellipse 

SLR  Y-coordinate  of  center  of  circle;  or  eccentricity  parameter  B 
in  equation  of  ellipse 
XCEN:  Xc  in  equation  of  ellipse 
YCEN:  Yc  in  equation  of  ellipse 

XR:  radius  of  circle 

Equation  of  ellipse: 

X-X-2  Y-Y2 

- -P-  +  —z^s - 1 


An  arc  which  begins  parallel  to  the  positive  X-axis  and  ends  parallel  to  the  positive  Y-axis 
has  a  starting  angle  of  3.14159  and  a  terminal  angle  of  1.57080. 

LT  <  0  signals  TSETUP  to  generate  a  dummy  tracer  particle  with  coordinates  (-1000,0)  to 
mark  the  end  of  a  tracer  string  for  a  package  or  subpackage. 

The  segments  must  form  a  continuous  boundary  except  where  they  intersect  the  grid 
boundaries.  Note  also  that  each  interface  is  described  by  two  identical  sets  of  tracer  particles, 
so  the  number  of  tracers  specified  along  an  interface  must  be  the  same  for  each  package 
although  generated  in  opposite  directions. 

Generally,  there  should  be  at  least  two  tracers  per  cell  edge.  Tn  regions  of  great  distortion 
tracer  particles  may  become  sparse;  however,  activating  ADDTCR  (see  Z-variable  NADD) 
helps  to  ensure  that  the  tracer  particles  do  not  become  too  sparse.  The  Z-variable  NIFMX 
limits  the  number  of  tracers  for  a  given  package  boundary.  NTPMX  should  equal  the  sum  of 
the  tracer  particle  arrays  TX  and  TY. 
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G.  Sipline  Endpoints 


Tf  sliplines  are  being  generated  (NOSLTP— 0)  SETUP  reads  data  defining  slipline  endpoints 
into  NAMELIST  SLIP.  Each  material  is  specified  as  a  master,  slave,  or  neither.  All  of  the 
SLIP  variables  can  be  read  in  a  single  SLIP  stream.  These  variables  are  defined  below. 


MASTRD 

NSLAVD 

NBGMD 

NBGSD 

NENDMD 

NENDSD 


—  N  if  package  (of  material  N)  is  a  master 

—  0  if  package  is  a  slave 

—  N  if  package  is  a  slave 

—  0  if  package  is  a  master 

=  index  of  first  tracer  particle  of  the  slipline  for  master 

—  0  for  slave 

—  0  for  master 

*=  index  of  first  tracer  particle  of  the  slipline  for  slave 
~  index  of  last  tracer  particle  of  the  slipline  for  master 

—  0  for  slave 
=  0  for  master 

=  index  of  last  tracer  particle  of  the  slipline  for  slave 


Sliplines  can  be  redefined  on  a  restart,  with  NOSLTP =0;  otherwise,  set  NOSLIP— 1  on  the 
restart. 


H.  High  Explosive  Detonation  Points 

SETUP  reads  data  defining  primary  and  secondary  initiation  points  into  NAMELIST  DETN. 
These  are  defined  below. 


IDET 


NDETS 

XDETS 

YDETS 

TIME 

XMIN 

YMTN 

XMAX 

YMAX 


=  1  for  primary  initiation  point 

—  2  for  secondary  initiation  point 
=  0  for  end  of  initiation  point  data 

—  material  package  number  of  high  explosive 

—  X-coordinate  of  initiation  point 

—  Y-coordinate  of  initiation  point 

—  delay  time  if  any 

=  left  X-coordinate  of  search  area 
=  lower  Y-coordinate  of  search  area 
=  right  X-coordinate  of  search  area 

—  upper  Y-coordinate  of  search  area 


The  first  five  of  these  variables  are  read  into  DETN  once  for  each  initiation  point.  TIME 
corresponds  to  the  delay  time  if  any  for  primary  initiation  points.  The  secondary  points  must 
be  read  in  the  order  in  which  detonation  will  occur.  The  remaining  variables  are  read  into 
DETN  once  for  each  detonation  point  to  define  a  search  area  for  the  code  to  determine 
detonation  times. 

A  card;  SDETN  IDET— 0$  must  be  read  in  problems  without  explosives. 
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I.  Restarts 


Restarts  require  only  the  NAMELIST’S  OPTNS,  START,  and  RUN  unless  NOSLIP— 0  in 
which  case  NAMELTST  SLIP  is  also  read. 

On  restarts  NAMELIST  RUN  allows  the  user  to  redefine  any  of  the  Z-block  -variables. 
However,  on  runs  which  terminated  on  TIMMAX  it  may  not  be  neccessary  to  change  any 
variable  including  TSTOP  or  ICSTOP,  depending  on  whether  the  specified  values  were  reached 
during  the  previous  run 


6,  PURE  CELLS  OVEREMPTYING 

Subroutine  TPHASE  determines  whether  a  pure  cell  has  overemptied  by  checking  the  sign 
of  the  variable  AMX(K),  where  AMX(K)  is  the  total  mass  in  the  Kth  cell.  Previously,  a 
negative  value  of  AMX  would  trigger  an  error  exit.  The  current  coding  conserves  mass  locally 
while  setting  the  mass  of  the  Kth  cell  to  zero.  Thus  the  code  is  allowed  to  continue  in  the  case 
where  the  overemptying  is  caused  by  small  numerical  error  during  a  stretching  condition. 

Tn  anticipation  of  some  other  problem  causing  the  overemptying,  a  diagnostic  message  is 
printed  at  the  time  of  the  recovery. 

The  modification  (TDENT  OVERMT)  is  shown  as  part  of  the  UPDATE  input  stream  in 
Appendix  A. 


7.  PRESSURE  ITERATION  PROBLEMS  IN  OLD  VERSIONS  OF  HELP 

For  cells  containing  more  than  one  material,  the  old  HELP  code  determined  cell  pressure 
by  an  iteration  procedure.  The  densities  of  the  separate  materials  in  such  a  cell  were  varied, 
with  specific  internal  energy  held  constant,  until  the  partial  pressures  converged  to  a  common 
value.  The  procedure  was  divided  into  a  pre-iteration  calculation,  to  ensure  that  the  cell  was 
filled  exactly,  followed  by  the  iteration  calculation. 

For  the  mixed  cell  it  was  required  that  the  partial  pressures  satisfy 


Pi(Pi,Vj)-P(p,V) 


where  P(p,V)  is  the  cell  pressure. 

In  the  iteration  process,  the  Vj  were  varied  until  the  P;  were  within  a  specified  tolerance  of  a 
common  pressure. 

Unfortunately,  the  pressure  iteration  routine  often  failed  to  converge  in  cases  where  the 
component  material  pressures  were  initially  approximately  equal,  but  the  relative  cell  volumes 
were  quite  different. 


14 


One  reason  why  the  convergence  routine  failed  is  that  the  JWLn  EOS  has  a  pressure 
maximum  which  is  potentially  hazardous  to  such  a  pressure  iteration.  The  EOS  may  be  called 
during  the  iteration  process  with  a  pressure  above  this  maximum.  Such  an  occurrence  would 
cause  failure  to  converge.  To  avoid  this,  the  variable  RHOMAX  was  defined  in  a  DATA 
statement  in  subroutine  CDT  as  the  maximum  density  to  be  allowed  for  an  HE  during  the 
pressure  iteration.  (RHOMAX— 9  for  Comp  B.)  However,  this  procedure  itself  often  caused 
the  convergence  routine  to  fail. 

Temporary  coding  allowed  recovery  by  restoration  of  the  initial  values  of  masses  and 
pressures  in  the  component  materials.  The  values  for  the  material  which  occupied  the  largest 
volume  in  the  cell  were  then  taken  as  the  mixed  cell  values.  Consequently,  a  warning  message 
was  primed  indicating  pressure  iteration  failure  and  recovery.  Material  pressures  for  each 
material,  before  and  after  the  iteration  attempt,  were  also  printed.  This  coding  allowed  the 
study  of  mixed-cell  conditions  which  led  to  pressure  iteration  failure. 

Table  1  shows  output  from  an  actual  calculation  in  which  the  pressure-iteration  failure 
recovery  was  invoked.  Pressures,  before  and  after  iteration,  are  shown  for  three  consecutive 
calculational  cycles.  The  pressures  for  each  material  are  approximately  equal  before  iteration, 
and  pressures  and  densities  for  one  or  both  materials  are  varied,  with  the  pressures  diverging 
during  the  iteration  process. 


TABLE  1 

Pressure  Iteration  Failure  Recovery  Data 


CY 

Pi 

Pi 

Pi 

P2 

Pi  IPi 

1055 

115.9 

2.25E+14 

8.66 

2.24E+14 

13.38 

121.9 

2.49E+14 

8.66 

2.24E+14 

1056 

115.8 

2.25E+14 

8.66 

2.24E+14 

13.37 

127.9 

2.75E+14 

9.56 

2.47E+14 

1057 

115.8 

2.25E+14 

8.66 

2.23E+14 

13.37 

133.8 

3.01E+14 

8.66 

2.23E+14 

77.  E.  Lee,  M.  Finger.  cmdW.  Collins,  'ML  Equation  of  State 
Coefficients  for  High  Explosive  s",  Lawrence  Livermore  Laboratory, 
UCID-161S9.  Jan  1973. 
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Table  2  shows  the  masses  and  mass  ratios,  at  5  cycle  intervals,  for  the  two  materials. 
During  the  period  shown,  material  1  is  passing  out  of  the  zone  and  material  2  is  occupying  an 
increasing  portion  of  the  total  cell  volume. 


TABLE  2 

Masses  and  Mass  Ratios  in  Mixed  Cells 


CY 

mj 

m2 

m,/m2 

1055 

6.67E-04 

8.59E-02 

7.76E-03 

1060 

5.42E-04 

1.06Er01 

5.11E-03 

1065 

4.04E-04 

1.25E-01 

3.23E-03 

1070 

2.68E-04 

1.42E-01 

1.89E-03 

1075 

1.10E-04 

1.51E-01 

7.28E-04 

1080 

6.50E-05 

1.54E-01 

4.22E-04 

1085 

3.34E-05 

1.55Er01 

2.15E-04 

1090 

3.86E-06 

1.55E-01 

2.49E-05 

Clearly,  the  initial  pressure  of  material  2  at  the  beginning  of  the  process  is,  to  a  good 
approximation,  the  correct  mixed-cell  pressure.  First,  the  two  material  pressures  are 
approximately  equal  and  second,  material  2  occupies  a  dominant  portion  of  the  cell  volume.  In 
such  a  case  the  iteration  process  itself  is  unnecessary  to  perform.  Furthermore,  it  is  apparent 
from  the  extremely  high  densities  of  material  1  that  the  iteration  construct  had  been  driving  the 
material  densities  to  unphysical  values,  well  before  the  failure  to  converge  had  indicated  a 
problem. 

Further  studies  with  3-material  mixed  cells  led  to  the  conclusion  that  a  mass-weighted 
average  pressure  for  mixed  cells  was  the  more  efficient  approach.  The  coding  which  invokes 
the  new  average  pressure  calculation  is  discussed  in  the  next  section. 
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8.  MASS-WEIGHTED  AVERAGE  PRESSURE  FOR  MIXED  CELLS 


Before  calculating  an  average  pressure  for  a  mixed  cell,  a  cell  filling  calculation  is  necessary 
because,  in  general,  the  masses  at  the  beginning  of  the  calculation,  do  not  exactly  fill  the  cell 
since  cell  masses  are  changed  in  the  TPHASE  portion  of  the  previous  calculational  cycle. 


The  code  computes  component  material  pressures  P(  from  component  densities  p|  while 
holding  the  specific  internal  energies  E)  constant  The  pj  are  then  varied  by  one  percent  to 
compute  new  Pj  and  constant  energy  compressibilities: 


Q2 


APi_ 

Ap,‘ 


New  volumes  are  calculated  using 


AVj 


The  pj  are  then  equated  to  1/Vj  where  the  total  cell  volume  is 


V-ZVi 


After  the  cell  filling,  the  mass-weighted  average  cell  pressures  are  computed  directly  from 


where 


and 


Ps  =  ZPimi 
Ws  =  5>i 


where  the  mj  are  component  masses. 
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9.  FTN  COMPILER  ERROR 


The  current  FTN  OPT1  compiler  at  the  BRL  computer  site  does  not  treat  certain 
NAMELIST  variable  addresses  correctly.  Apparently,  certain  three  letter  variable  addresses  are 
not  known  to  have  the  same  address  by  the  NAMELIST  package  and  the  rest  of  the  code.  In 
particular,  values  of  the  parameters  RMU(I),  the  rigidity  moduli  of  the  materials,  the  UUR(l), 
initial  radial  (or  x)  velocities,  and  WA(I),  initial  axial  (or  y)  velocities,  read  as  non-zero  into 
the  location  reserved  for  them  by  NAMELIST,  remained  zero  to  the  rest  of  the  code  because 
the  variables  of  the  same  name  known  to  the  rest  of  the  code  had  different  addresses. 

For  the  duration  of  this  compiler  residing  on  the  BRL  Cyber  system,  corrective  coding 
(IDENT  MAY80FTNE),  shown  in  Appendix  A,  is  neccessary  for  avoiding  this  problem. 


10.  BRLEOS  CODE  SYSTEM 

In  old  versions  of  HELP  the  subroutine  EQST  performed  thermodynamic  calculations  using 
the  Tillotson  equation  of  state  for  metals  and  composites,  the  JWL  equation  of  state  for  explo¬ 
sives,  and  an  ideal  gas  equation  of  state  for  gases. 

In  the  new  version  of  the  code,  subroutine  EQST  is  used  to  initialize  certain  EOS  data  and 
call  subroutine  BRLEOS.  BRLEOS  includes  all  of  the  EOS  coding  of  the  old  EQST  plus  calls  to 
the  BRLGRAY  equation  of  state  and  the  BRLPUFF  equation  of  state. 

Applications  of  the  EOS’s  are  listed  in  Table  3  and  material  index  numbers  are  listed  in 
Table  4. 


TABLE  3 

Equation  of  State  Applications 

EOS 

Application 

Type 

BRLGRAY 

Metals/ Alloys 

Mixed-Phase 

Solid-Liquid-Vapor 

Tillotson 

Metals/Composites 

Joined  Compressed-Exp anded 

Gamma-Law 

Ideal  Gas 

P={y-l]pE 

JWL 

Explosives 

P(E,p) 

BRLPUFF 

Metals/Phenolics 

Joined  Compressed-Expanded 
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TABLE  4 

Equation-of-State  Material  Indices 

Index 

Material 

EOS 

1 

Aluminum 

2 

Iron 

3 

Copper 

4 

Tantalum 

5 

Uranium 

6 

Magnesium 

7 

Berylium 

8 

Titanium 

9 

Tin 

BRLGRAY 

10 

Stainless 

11 

Niobium 

12 

Nickel 

13 

Molybdenum 

14 

Lead 

15 

Thorium 

16 

Tungsten 

17 

Gold 

18 

RHA 

101 

Tungsten 

102 

Copper 

103 

Iron 

104 

Aluminum 

105 

Berylium 

106 

Titanium 

107 

Nickel 

108 

Molybdenum 

109 

Thorium 

110 

Lead 

Tillotson 

111 

Ch2 

112 

Granite 

113 

Andesite 

114 

Wet  Tuff 

115 

Dry  Tuff 

116 

Oil  Shale 

117 

Dolomite 

118 

Limestone 

119 

Halite 

120 

Ideal  Gas 

Gamma-Law 

121 

Comp  B 

122 

TNT 

123 

Octol 

124 

PBX9404 

JWL 

125 

LX-14-0 

126 

Hivelite 

227 

Aluminum 

228 

Plexiglas 

BRLPUFF 

229 

Fused  Quartz 

230 

Quartz  Phenolic 

19 


The  EOS  constants  for  materials  112  to  119  were  given  by  Evans  and  Harlow12  and  it  should  be 
noted  that  the  constants  in  the  Tillotson  EOS  for  the  materials  Granite  and  Andesite  are  valid 
only  for  pressures  less  than  150  kbar. 

11.  SAMPLE  CALCULATION 

A.  Sample  Input 

The  jobstream  corresponding  to  Sample  Problem  2  of  Reference  1  is  shown  in  Appendix  R 
As  the  first  input  card  indicates,  this  calculation  is  of  a  copper,  hemispherical  liner  loaded  with 
unconfined  Comp  B. 

After  the  edit  control  parameters  and  Z-variables  are  read,  the  new  NAMELIST  input 
begins.  The  grid  that  is  specified  has  60  by  200  square  cells  of  .05cm  dimension.  The  next 
card  specifies  that  the  cell  in  the  bottom  left  corner  of  the  grid  is  a  pure  void  cell.  The  rest  of 
the  MTRL  input  specifies  material  1  as  copper  and  material  2  as  Corap  R 

The  GEOM  input  begins  with  the  specification  of  an  arc  of  a  circle  (LT“4)  bounding 
material  1  with  300  tracers.  The  starting  angle  is  PI  and  the  terminal  angle  is  PT/2.  The  center 
of  the  circle  is  at  (0.0,6.35)  with  a  radius  of  1.905.  The  next  segment  is  a  horizontal  line  going 
from  (1.905,6.35)  to  (1.7907,6.35)  with  50  tracers.  The  description  of  the  liner  is  completed 
with  an  arc  (LT=-4)  of  radius  1.7907  corresponding  to  the  inner  radius  of  the  liner. 

The  remainder  of  the  GEOM  input  bounds  the  Comp  B  charge  (MPN=2)  and  the  void 
(MPN— 3). 

In  this  problem  a  slipline  is  placed  along  the  interface  between  the  liner  and  the  explosive. 
The  SLIP  input  is  straightforward.  The  indexing  begins  with  the  master  material  on  the  left 
starting  with  the  first  tracer  particle  of  the  master  material,  counting  up  to  the  last  master  tracer 
particle  on  the  interface  (300),  counting  back  along  the  interface  with  the  slave  material  on  the 
left  starting  with  the  first  slave  tracer  (301)  and  ending  with  the  last  slave  tracer  on  the 
interface  (600). 

For  this  explosive  configuration,  both  a  primary  and  secondary  detonation  point  are 
required.  They  are  placed  at  (0.0,2.54)  and  (1.82,5.62),  respectively.  Zero  delay  time  is 
specified  in  both  cases.  The  search  areas  specified,  for  the  code  to  calculate  cell  detonation 
times,  have  lower  left  and  upper  right  corners:  (0.0,2.54)  and  (2.159,6.35),  and  (0.0,5.59)  and 
(2.159,6.35),  respectively.  Less  overlap  would  be  desirable  but  the  computing  time  involved  is 
not  significant.  What  is  important  is  the  placement  of  the  secondary  initiation  point  which 
allows  it  to  "shine"  on  the  portions  of  the  explosive  not  covered  by  the  primary  point 


12.  M.  W.  Evans  and  F.H.  Harlow.  "The  Particle-ln-CeH  Method Jbr 
Hydrodynamics  Calculation^,  Los  Alamos  Scientific  Laboratory 
Report  LA-2139,  Nov  19 57. 
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B.  Sample  Output 


Appendix  C  contains  the  output  for  Sample  Problem  2  run  on  BRLHELP,  cyde  23.  The 
first  page  shown  is  part  of  the  dayfile.  A  copy  of  a  user-library,  LIBCMC,  containing  the 
BRLEOS  codes  is  attached,  along  with  the  main  program  library  BRLHELP,  and  the 
MIXEDPHASEBOUNDARY  data  tape  required  by  BRLGRAY.  This  was  a  short  run,  using 
867.305  JOB  seconds  (corresponding  to  TTMMAX—900).  The  output  was  stored,  for  later  use 
on  a  restart,  on  a  permanent  disk  file  called  CUHEMT. 

The  first  page  of  actual  code  output  lists  the  edit  parameters.  The  next  page  shows  the 
problem  identification  (read  on  the  header  card)  and  the  geometry  package  including  definition 
of  the  slide-line  endpoints.  This  is  followed  by  a  description  of  the  explosive  package  on  two 
consecutive  pages.  The  Z-variables  are  printed  next,  followed  by  a  general  description  of  the 
problem  grid. 

The  material  in  package  1,  G  Copper,  is  copper  described  by  the  BRLGRAY  EOS.  The 
material  in  package  2,  J  Comp  B,  is  Comp  B  described  by  the  JWL  EOS.  Next  comes  the 
strength  constants  and  a  description  of  the  grid.  This  takes  several  pages  of  output  concluding 
with  the  cell-coordinates  of  the  tracer  panicles  for  each  material  package. 

The  compression  map  at  cycle  0  shows  the  compression  values  of  both  the  metal  liner  and 
the  explosive  (mixed-cells  are  indicated  by  asterisks)  to  have  the  symbolic  value  21  Z  in  this 
case  has  the  value  of  1.0  where  in  general  the  value  of  Z  is  given  as  a  range  of  values,  for 
example  from  .980  to  1.026.  The  display  of  mixed  and  pure  cells,  shown  next,  shows  the  metal 
liner  as  material  1  and  the  explosive  as  material  2,  with  mixed  cells  denoted  symbolically  with 
the  letter  M.  Note  that  the  void  is  considered  as  a  material  for  computational  purposes. 

The  next  page  shown  is  the  end  of  the  cycle  0  printout  and  the  beginning  of  the  regular 
cyclic  printout  This  continues  until  it  is  time  for  a  tape  dump  and  printout  at  a  specified  time 
or  cycle,  or  at  the  end  of  the  calculation.  In  this  case  the  next  tape  dump  was  at  the  specified 
time  of  2.0  /ns  which  occurred  at  cycle  69. 

The  next  page  shown  is  the  pressure  map  at  2.0  /xs.  The  detonation  wave  has  not  quite 
reached  the  metal  liner,  and  the  BRLGRAY  EOS  has  not  yet  been  called  for  a  pressure  and 
temperature  calculatioa  When  the  shock  wave  does  hit  the  liner  and  the  EOS  is  called  for  the 
first  time,  it  prints  out  a  set  of  input  (ESTCON)  and  calculated  (EGONS)  equation-of-state 
constants.  This  occurs  during  the  regular  cyclic  printout 

The  last  page  of  output  shown  is  the  temperature  map.  Temperatures  in  the  liner  range 
symbolically  from  E  to  H,  indicating  temperatures  from  416  to  832  degrees  Kelvia  The  Q’s 
and  the  R  occurring  along  the  explosive/metal  interface  are  transient  values  which  have  little 
effect  on  the  calculatioa 
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12.  SUMMARY 


The  HELP  code  has  been  shown  to  be  an  excellent  tool  for  shaped-charge  liner  optimization 
studies;  however  there  have  been  some  nagging  problems  for  its  users.  Two  of  these  have 
been  addressed.  These  are  the  fatal  errors  caused  by.  (1)  over-emptying  of  pure  cells  and  (2) 
the  failure  of  the  pressure  iteration  routine  to  achieve  convergence  for  certain  mixed-cell 
situations.  A  cycle  23  version  of  BRLHELP  is  now  available  which  alleviates  these  problems. 
The  overemptying  of  pure  cells  has  been  converted  to  a  non-fatal  condition  and  pressure 
iteration  failure  for  mixed  cells  has  been  eliminated  by  employment  of  a  mass-weighted  average 
pressure  scheme. 

A  call  to  the  BRLEOS  code  system  has  been  installed  in  place  of  the  old  EQST  subroutine  to 
allow  access  to  a  variety  of  equations  of  state  including  BRLGRAY  which  provides  temperature 
as  an  output  variable.  A  problem  with  the  Fortran  compiler  on  the  BRL  Cyber  176  has  been 
circumvented. 

Input  to  the  BRLHELP  code  has  been  changed  to  the  NAMELIST  format.  Interpretation 
of  input  decks  is  thereby  less  cryptic  and  modifications  to  old  decks  less  prone  to  error. 

13.  RECOMMENDATIONS 

BRLHELP  should  be  considered  for  more  extensive  use,  at  the  BRL,  in  predictive  studies. 
Its  new  capability,  via  BRLGRAY,  to  produce  temperature  as  an  output  variable  should  prove 
most  useful  to  researchers  trying  to  determine  the  state  of  shocked  materials.  Its  new 
robustness  should  prove  the  code  to  be  especially  useful  in  prefiring  analyses,  where  use  of  the 
code  can  save  considerable  expense  by  avoiding  trial  and  error  experimentation. 


11-  E.  Lee ,  M .  Finger,  andW .  Collins ,  *  JW  L  Equation  of  State 
Coefficients  for  High  Explosive  s*,  Lawrence  Livermore  Laboratory, 
UCID-161K9,  Jan  1973 . 

12.  M.  W.  Evans  and  F.H.  Harlow ,  "The  Partlcle-fn-Cell  Method  for 
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APPENDIX  A:  UPDATE  LISTINGS 
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I£TRA,STWFZ,T200.  UPDATE  TO  BRIHELP;  CY=17  TO  CY=18 
ACCOUNT (MD440) 

REQUEST (B ,  *PF) 

GETPF (A,BRLHELP,CY=17, ID=CMCHYDRO) 

UPDATE (F,P=A,N=B) 

FTN (I=COMPILE,L=0,PL=55000,ICM=I) 

CATALOG (B,BRLHELP ,CY=8 , ID=CMCHYDRO) 

*EOR 

*IDENT  OVERMT 

*D  TPHASE . 687 , TPHASE .689 

C  ***  IF  PURE  CELL  OVEREMPTIED  CONSERVE  MASS  LOCALLY  *** 

C  ***  AND  ZERO  OUT  OVEREMPTIED  CELL.  *** 

1000  IF  (AMX (K) )  1010,1070,1015 
1010  CONTINUE 

WRITE (6,1310)  I,J 
KKA=K-IMAX 

WRITE (6,1350)  I,J,K,MFIAG(K)  ,AMX(K)  ,AIX(K)  ,U(K)  ,V(K)  ,DETIM (K)  , 
1  MFIAG(KKA)  ,AMX  (KKA)  ,AIX(KKA)  ,U(KKA)  ,V(KKA)  ,DETIM  (KKA)  ,ETH 
1350  FORMAT (3I5,/I5,5E12.5,/I5,6E12.5) 

AMX (KKA) =AMX (KKA) -AMK  (K) 

AMX (K)=0. 

GO  TO  1060 
*IDENT  LACITERA 
*1  CDT.12 

DIMENSION  DNSAV ( 4 ) ,PRSAV(4) ,XMSAV(4) 

*1  CDT.71 

DNSAV(L)=0. 

PRSAV(L)=0. 

XMSAV  (L)  =0 . 

VMAX=0 . 

PMXVL=0. 

*1  CDT.82 

DNSAV(L)=RHO(L,M) 

XMSAV  (L)=XMASS  (L,M) 

*1  CET.129 

PRSAV (L) =PRS (2,L) 

IF (VMAX.GT. VOL (L) )  GO  TO  120 
VMAX=VOL (L) 

PMXVL=PRS (2 ,  L) 

*D  CDT.246 

NK=260 

WRITE (6,560)  NK 
*D  CDT. 249,CDT. 250 
GO  TO  314 
*1  CDT. 283 
301  CONTINUE 
*D  CDT. 301 
NK=310 

WRITE (6,560)  NK 
*D  CET. 305, CDT. 306 
314  DC  315  L=1,NMAT 
XMASS  (L,M)=XMSAV(L) 

DNS  (1 ,  L)  =DNSAV  (L) 

PRS  (1 ,  L)  =PRSAV  (L) 
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LCTRA,STMFZ,T200.  UPDATE  TO  BRIHELP;  CY=17  TO  CY«=18 
ACCOUNT (MD440) 

REQUEST (B,  *PF) 

GETPF  (A,BRLHELP,CY=17,  ID=CMCHYDRO) 

UPDATE (F,P=A,N=B) 

FIN  (I COMPILE ,  L=0 , PL=55000 , LCM=I ) 

CATALOG  (B ,  BRLHEIP  ,CY=8 ,  IDCMCHYDRO) 

*EOR 

*IDENT  OVERMT 
*D  TPHASE.687,TPHASE.689 

C  ***  IF  PURE  CELL  OVEREMPTIED  CONSERVE  MASS  LOCALLY  *** 

C  ***  AND  ZERO  OUT  OVEREMPTIED  CELL.  *** 

1000  IF  (AMX (K) )  1010,1070,1015 
1010  CONTINUE 

WRITE (6,1310)  I,J 
KKA=K-IMAX 

WRITE (6,1350)  I,J,K,MFIAG(K)  ,AMX(K)  ,AIX(K)  ,U(K)  ,V(K)  ,DETIM(K) 
1  MFIAG(KKA)  ,AMX (KKA)  ,AIX(KKA)  ,U(KKA)  ,V(KKA)  ,DETIM (KKA)  ,ETH 
1350  FCRMAT(3I5,/I5,5E12.5,/I5,6E12.5) 

AMX  (KKA)  =AMX  (KKA)  -AMK  (K) 

AMX  (K)  =0 . 

GO  TO  1060 
*IDENT  LACITERA 
*1  CDT.12 

DIMENSION  DNSAV(4) ,PRSAV(4) ,XMSAV(4) 

*1  CDT.71 

DNSAV (L) =0 . 

PRSAV(L)=0. 

XMSAV (L) =0 . 

VMAX=0. 

PMXVL=0. 

*1  CDT.82 

DNSAV  (L)=RHO(L,M) 

XMSAV (L)=XMASS  (L,M) 

*1  CDT.129 

PRSAV  (L)  =PRS  (2 ,  L) 

IF (VMAX.GT.VOL(L) )  GO  TO  120 
VMAX=VOL (L) 

PMXVL=PRS ( 2 , L) 

*D  CDT.246 

NK=260 

WRITE (6,560)  NK 
*D  CDT. 249 ,CDT. 250 
GO  TO  314 
*1  CDT. 283 
301  CONTINUE 
*D  CDT. 301 

NK=310 

WRITE (6,560)  NK 
*D  CDT. 305, CDT. 306 
314  DO  315  L=1,NMAT 
XMASS  (L,M)  =XMSAV  (L) 

INS  (1 ,  L)  =DNSAV  (L) 

PRS  (1 ,  L)  =PRSAV  (L) 
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IF (PRS (1,L) .GT.PMIN)  PSUM=1 
315  PAV=PMXVL 
*CCPY  CDT,CDT.303,CDT.304 
GO  TO  301 
*D  CDT.439 

560  FORMAT (IHI , *PRESSURE  ITERATION  FAILURE,  RECOVERED  AT  CUT 

1  LINE  NUMBER*, 15) 


ICTRA ,  STWFZ ,  T200 .  UPDATE  TO  BRLHELP;  CY=18  TO  CY=I9 
ACCOUNT (MD440) 

REQUEST (B,  *PF) 

GETPF (A , BRLHELP , ID=CMCFYDRO) 

UPDATE (F,P=A,N=B) 

PIN  (I=COMPILE , L=0 , PL=55000 , LCM=I ) 

CATALOG (B, BRLHELP, ID=CMCHYDRO,RP=999) 

*EOR 

*IDENT  MAY80FTNE 
*D  GEOM.7 

I  CZERO,  STKI ,  STK2,  STEZ,  TEMP,  AMDM 

*1  GECM.IO 

DO  100  IRMU=1 , 4 
100  RMU  (IRMU)  =TEMP  (IRMU) 

C  ***  THIS  CODING  NECESSARY  TO  AVOID  CPTl  COMPILER  ERROR  *** 


LCTRA,STMFZ,T200.  UPDATE  TO  PRIHELP;  CY=19  TO  CY=20 
ACCOUNT (MD440) 

REQUEST (B,  *PF) 

GETPF  (A,  BRIHELP ,  ID=€MCHYDRO ) 

UPDATE  (F,P=A,N=B) 

PIN  (I=CCMPILE ,  L=0 ,  PL=55000 ,  ICM=I ) 

CATALOG  (B, BRIHELP,  ID=CKCHYDRC,RP=999) 

*EOR 

*IDENT  CELLFILL 
*YANK  IACITERA 
*D  CDT.78 
*D  CDT.95,CDT.162 
C 

C  ***  COMPUTE  TWO  DENSITIES  FOR 
C  ***  EACH  MATERIAL  IN  CELL  USING 
C  ***  (1)  DENSITIES  FRCM  LAST  CYCLE 

C  ***  (2)  THOSE  DENSITIES  INCREASED  BY  ONE  PERCENT 

C 

*COPY  CDT,CET.99,CDT.129 
*COPY  CDT,CDT.130,CDT.141 

IF  (VL(L).CT.O)  GO  TO  170 
VL(L)=0.0 
XMASS (L,M)=0.0 
170  CONTINUE 
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*IDENT  PCALC 

*D  CUT. 163, 316 

*COPY  CDT,CDT.163,CDT,176 

*COPY  CDT,CDT. 184 ,CDT. 224 

C 

C  ***  MASS  WEIGHTED  AVERAGE  PRESSURE 
P(K)=PSUM/WSUM 

C 

*COPY  CDT,CDT.308,CDT.316 
*IDENT  LCTSEP80 
*D  INPUT. 75, INPUT. 76 
G 


CALL  SETUP 

C  ***  ONLY  TO  READ  SLIP  LINE  DATA  AND  RETURN  *** 
C 

*B  SETUP. 13 
C 

C  ***  ip  RESTART  READ  SLIP  DATA  AND  RETURN  *** 

IF (PK (3) .LT.O)  READ (5, SLIP) 

IF (PK (3) .LT.O)  RETURN 
C 

*D  GEOM.6,MAY80FTNE.l 
C 


NAMELIST/MTRL/  MFIAG,  MAT, 

1  CZERO,  STK1 ,  STK2, 

*D  MAY80FTNE .  2  ,MAY80FTNE .  4 
DO  100  1=1,4 


RHOIN,  SSIEN,  UURX,  WAY, 
STEZ ,  RMUI ,  AMDM 


UUR (I)=UURX (I) 
WA  (I )  =WAY  (I ) 
100  RMU (I)=RMUI (I) 


C  ***  THIS  CODING  NECESSARY  TO  AVOID  OPTl  COMPILER  ERROR* 

C  ***  DETECTED  BY  JOE  IACETERA  DURING  PERIOD  MAY-SEP  198( 
*1  SETUP. 11 


*** 

*** 


DIMENSION  UURX(4) ,  WAY (4),  RMUI  (4) 


ICTRA,STMFZ,T100.  UPDATE  TO  BRIHELP;  CY=20  TO  CY=21 
ACCOUNT (MD440) 

REQUEST (NEWPL,  *PF) 

GETPF  (OIDPL  ,ERLHELP ,  ID=CMCHYDRO) 

UPDATE (F,N) 

CATALOG  (NEWPL, BRIHELP,ID=CMCHYDRD) 

EXIT. 

*EOR 

*IDENT  EQSTEOS 
*D  BRL76.6 

1  TAPE10 , TAPE1 ) 

*1  EQST.6 


*1 


COMMCN/EQSTEOS/XAMDM (4 )  ,XEMIN,XGAMMA,XROEPS,XMFK 
HELPCOM. 71  ' 

I£VEL  2,  PLWP , PLWC , PVOL , ALPSV 
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*D  EQST.11,EQST.234 


XAMDM  (MFK)  «=AMDM  (MFK) 

XEWIN=EMIN 

XGAMMA=GAMMA 

XROEPS=ROEPS 

XMFK=MFK 

ENERGY=AMAX1 (ENERGY ,  0 . ) 

CALL  ERLEOS (ENERGY, RHOW,PRESUR,TWPTR,CS, IV,  N) 

*D  MEMEXPND.l 

COMMCN/MXCELL/  RHOZ (130) ,CNAUT(130) ,MAT(130) , SAMMY (4,800) , 
*D  INPUT. 8 

DIMENSION  MNAME (130) 

*D  HELPEQ.18,HELPEQ.36 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


MATERIAL  PROPERTY  DATA  STATEMENTS 
MATERIALS  I  -  18  ARE  ERIGRAY  EOS  METALS 
MATERIALS  101-119  ARE  TILLCTSON  MATERIALS 
MATERIAL  120  IS  IDEAL  GAS 

NOTE  -  IDEAL  GAS  REFERENCE  DENSITY  RHOZ  (120)  IS  THE 
NORMAL  DENSITY  OF  AIR.  THIS  VALUE  MUST  BE 
CHANGED  FOR  OTHER  GASES. 

EXPLOSIVE  121  IS  CCMPB 
EXPLOSIVE  122  IS  TOT 
EXPLOSIVE  123  IS  OCTOL 
EXPLOSIVE  124  IS  PBX  9404 
EXPLOSIVE  125  IS  LX-14-0 
EXPLOSIVE  126  IS  HIVELITE 


DATA  (RHOZ(K) ,K=1,18)  / 


1 

2.70, 

7.87, 

8.94, 

16.65, 

19.08, 

2 

1.75, 

1.85, 

4.51, 

7.28, 

7.91, 

3 

8.58, 

8.86, 

10.22, 

11.36, 

11.68, 

4 

19.17, 

19.24, 

7.84/ 

DATA 

(RHOZ (K) 

,K=101 ,119) 

/ 

1 

19.17, 

8.94, 

7.87, 

2.70, 

1.85, 

2 

4.51, 

8.86, 

10.22, 

11.68, 

11.36, 

3 

.90, 

2.70, 

2.70, 

1.97, 

1.70, 

4 

2.30, 

2.80, 

2.70, 

2.20/ 

DATA 

RHOZ (120)/.001293/ 

DATA 

(RHOZ (K) 

,K=121,126) 

/ 

1 

1.717, 

1.630, 

1.821, 

1.840, 

1.835 
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DATA  (CNAUT(K) 

1  4.01E5, 

2  4.78E5, 

3  2.89E5, 

4  3.49E5, 


K=101,119)/ 
3.95E5,  4.03E5, 

4.63E5,  5.15E5, 

2.58E5,  2.58E5, 

5.51E5,  3.85E5, 


5.27E5,  8.06E5, 
2.13E5,  2.03E5, 
2.24E5,  1.63E5, 
3.37E5/ 


*D  CDT.15,CET.22 

GAMGAS  (120)  =  GAMMA 

DATA  (GAMGAS (K),  K=121,126)/2.706,2.727,2.830,2.850,2.841,1.24/ 
DATA  (RHCMAX(K) ,  K=121,126)/9.0,  8.5,  9.5,  11.0,  9.5,  9.5/ 

C 


*D  DETIME. 7, DETIME. 18 
C 

C  EXPLOSIVE  121  IS  CCMPB 

C  EXPLOSIVE  122  IS  TNT 

C  EXPLOSIVE  123  IS  OCTOL 

C  EXPLOSIVE  124  IS  PBX  9404 

C  EXPLOSIVE  125  IS  LX-14-0 

C  EXPLOSIVE  126  IS  BIVELITE 

C 

C  DETVEL  IS  THE  DETONATION  VELOCITY 
C 


DATA  DETVEL (121)  /7.98E5/ 
DATA  DETVEL (122)  /6.93E5/ 
DATA  DETVEL (123)  /8.48E5/ 
DATA  DETVEL (124)  /8.80E5/ 
DATA  DETVEL (125)  /8.83E5/ 
DATA  DETVEL (126)  /.381E5/ 


*D  ADDENG . & , ADDENG . 1 2 

DATA  (DENGY(K) ,K=121,126)  / 

2  4.95E+10,  3.68E+10,  5.27E+10,  5.54E+10,  1.02E+10,  4.20E+10/ 

C  EXPLOSIVE  121  IS  CCMPB 

C  EXPLOSIVE  122  IS  TNT 

C  EXPLOSIVE  123  IS  OCTOL 

C  EXPLOSIVE  124  IS  PBX  9404 

C  EXPLOSIVE  125  IS  LX-14-0 

C  EXPLOSIVE  126  IS  HIVELITE 

*D  INPUT.  10, INPUT.  16 

DATA  (MNAME(K) ,K=1,18)/ 

1  10HG  ALUMINUM, 10HG  IRON,10HG  COPPER, 10HG  TANTALUM, 

2  10HG  URANIUM, 10HGMAGNESIUM,10HG  BEFYLIUM , 10HG  TITANIUM, 

3  10BG  TIN,10HG  S-STEEL,10HG  NICBIUM,10HG  NICKEL, 

4  10HG  MOLY,10HG  LEAD,10HG  THORIUM, 10HG  TUNGSTEN, 

5  10HG  GOLD,10HG  RHA/ 


DATA  (MNAME(K) ,K=101,124)/ 

10HT  TUNGSTEN, 10HT  COPPER, 10HT 
10HT  BERYLIUM , 10HT  TITANIUM, 10HT 


1 

2 

3 

4 

5 

6 
7 


10HT  THORIUM, 10HT 
10HT  POLYMERS, 10HT 
10HT  DRY  TUFF,10HT 
1 OHTLIMESTCNE , 10HT 
10HJ  TNT,10HJ 


IRON,10HT 
NICKEL, 10HT 

LEAD, 

GRANITE, 10HT  ANDESITE, 10HT 
OIL,10HT  SHALE, 10HT 

HALITE, 10H  IDEAL  GAS,10BJ 
OCTOL, 10HJ  PBX, 1 OBJ 


ALUMINUM, 

MOLY, 

WET  TUFF, 
DOLOMITE, 
COMP  B, 
9404/ 
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DATA  MNAME (125)  /10HJ  LX-14-0/ 

DATA  MNAME (126)  /10HJ  HIVELITE/ 

*D  INPUT. 99, INPUT. 107 

C  PRINT  INITIAL  CONDITIONS 

WRITE (6,390) 

DO  60  L=1,NMAT 
MN=MAT  (L) 

WRITE  (6,400)  L,RHOZ  (MN)  ,RHOIN (L)  ,SSIEN  (L)  ,UUR  (L)  ,WA  (L) , MNAME  (MN) 
60  CONTINUE 
*D  ERL76.90 

400  FORMAT ( 1 5 , 2F13 . 3 , 9X , 1PE1 1 . 4 , 5X , 1PE11 . 4 , 5X , 1PE11 . 4 , A 10 ) 

*1  EDIT. 240 

305  CONTINUE 
*D  NOEDPRT.8 


LCTRA,STMFZ,T200.  UPDATE  TO  ERIEELP;  CY=21  TO  CY=22 
ACCOUNT (MD440) 

RECUEST (B ,  *PF) 

GETPF (TA , PLCMCEOS , ID=CMCeOS ) 

UPDATE (Q,P=TA,N=TEMPL) 

GETPF  (A ,  BRLHELP ,  33XMCHYDRO) 

UPDATE  (F,P=A,M=TEMPL,N=B) 

FIN  (I=COMPILE ,  L=0 , PL=55000 , LCM=I ) 

CATALOG  (B,ERIHELP,ID=CMCHYDRO,RP=999) 

*EOR 

♦COMPILE  BRIEOS 
*EOR 


*IDENT  ECSTEOSCR 
*D  EQSTEOS . 45 ,EQSTEOS . 47 

5  10HT  DRY  TUFF, lOlfr  OILSHAI£,10HT  DOLOMITE,  10HTLIMESTONE, 

6  10HT  HAUTE/ 

C 


C 


DATA  MNAME (120)  /10H  IDEAL  GAS/ 


DATA  MNAME (121)  /10HJ  COMP  B/ 

DATA  MNAME  (122)  /10HJ  TNT/ 

DATA  MNAME (123)  /10HJ  OCTOL/ 

DATA  MNAME (124)  /10HJ  PBX  9404/ 

♦IDENT  EOSDET 
*D  CDT.12 


DIMENSION  RHCMAX (130)  ,  GAMGAS(130) 

*D  DETIME. 5 

1DETVEL(130) 

*D  ADDENG . 4 , ADDENG. 6 

DIMENSION  DENGY (130) 

*D  DETIME. 132 

IF(N.GT.120.AND.N.U.200)  GO  TO  105 
*D  DETIME. 135 

100  IF(MAT(M) .LE.120.QR.MAT(M) .GE.200)  GO  TO  120 
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*D  DETIME. 212 

IF (MAT(MFK)  .GT. 120. AND. MAT (MFK)  . LT.200)  GO  TO  190 
*D  DETIME. 21 8 

IF (MAT (L) .GT. 120. AND.MAT(L) .LT.200)  GO  TO  190 
*D  EETIME. 235 

IF (MAT(L) .GT.120.AND.MAT(L) .LT.200)  GO  TO  173 
*D  DETIME. 24 5 

IF (MAT(L) .GT. 120. AND. MAT (L) .LT.200)  GO  TO  176 
*D  DETIME. 256 

IF (MAT(MFL)  .CT.120.AND. MAT (MFL)  .LT.200)  GO  TO  181 
*D  DETIME. 261 

IF (MAT (L).GT.120.AND.MAT(L). LT.200)  GO  TO  181 
*D  EETIME. 272 

IF (MAT(MFL)  .GT. 120.AND. MAT(MFL)  .LT.200)  GO  TO  185 
*D  EETIME. 276 

IF (MAT (L) .GT. 120. AND. MAT (L). LT.200)  GO  TO  185 
*D  DETIME. 365 

IF  (MAT(MFK)  .IE.  120. OR. MAT (MFK)  .(£.200)  GO  TO  170 
*D  DETIME. 393 

IF (MAT (L) ,I£.120.OR.MAT(L) .GE.200)  GO  TO  360 


LCTRA , STMFZ , T200 .  UPDATE  TO  BRIflELP;  CY=22  TO  CY=23 
ACCOUNT (MD440) 

PECUEST (B ,  *PF) 

GETPF (A,BRIHELP, ID=CMCHYDPO) 

UPDATE (F,P=A,N=B) 

FTN  (I=COMPILE ,  L=0 ,  PL=55000 ,  ICM=I ) 

CATALOG  (B ,  BRLHELP ,  ID=CMC  HYDRO ,  RP=9  9 9 ) 

*EOR 

*IDENT  GEOMCR 
*D  GECM.15 

30  WRITE (6,370)  X1S,YIS,X2F,Y2F 
*IDENT  IMPOCT80 
*1  HELPCOM. 83 

COMMCN/TGRAY/TGRAY,  TXG  (12000) ,  TRS(4) 

LEVEL  2,  TGRAY,  T)G,  TRS 

*D  E0ST.3,ECST.6 

C  ***  INITIALIZE  SCME  BRIGRAY  PARAMETERS  AND  CALL  BRIE  OS 
C 

*B  ECST.32 

TGRAY=0. 

*1  ECST.32 

TGRAY=TMPTR 
*B  CDT.23 

TGRAY=0 
*1  CDT.60 

TXG (K)=TGRAY 
*1  PCALC.12 

TRS  (L)  =TGRAY 
*1  PCALC.35 
TSUM=0. 

*1  PCALC.54 

TSUM=AMAXl  (TSUM,TRS  (L) ) 
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non 


*1  PCALC.58 
C 

C  ***  MAXIMUM  TEMPERATURE  FCR  ANY  METAL  IN  TOE  MIXED  CELL 
TXG(K)=TSUM 

*D  MAP. 7 

DIMENSION  WSMAX (6) ,  AI£(41),  XUM(41),  VALUE (41) 

*1  MAP. 71 

C  ***  TEMPERATURE 

DO  75  J=1,JDL 
DO  75  1=1 , IDL 
K= (J-l ) *IMAX+I+1 
M=IAES  (MFIAG(K)) 

N=MAT (M) 

IF (N.GE.100)  TXG (K) =0. 

IF (N.GE.100)  GO  TO  74 
IF  (TXG  (K)  .EQ.O)  TX3(K)=300. 

74  CONTINUE 

75  WSMAX  (6)=AMAX1  (WSMAX  (6)  ,ABS  (TXG  (K) ) ) 

PRINT  *,  "IMAX=" , WSMAX  (6 ) 

*1  MAP. 139 

GO  TO  230 

***  TEMPERATURE 


225  J-JDL 

IF  (WSMAX  (NPROP)  .LE.O)  GO  TO  460 
WRITE (6,535)  CYCLE, T 

226  DO  228  1=1,  IDL 
K= (J-l ) *IMAX+I+1 

228  PROP(I)=TXG (K) 

*D  MAP. 234 

GO  TO  (90,120,150,180,210,226)  NPROP 
*D  MAP. 242 

GO  TO  (470,110,140,170,200,225,470)  NPROP 
*1  MAP. 259 

535  FORMAT (1H1 , 4X , 15HTEMPERATORE  ,1X,6HCYCLE=,F8.1,4X,5HTIME=, 

*1PE11 . 5 , IX , 7HSECONDS//) 
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APPENDIX  B:  SAMPLE  INPUT 
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LCTIL,STMFZ,T1000,P0.  SAMPLE  PROBLEM  2:  CU  HEMISPHERE  LACETERA 
ACCOUNT (MD440) 

REQUEST  (TAPE7 ,  *PF) 

GETPF (LIBCMC , LIBCMC ,  ID=CMCEOS) 

GETPF (A ,  BRLHELP ,  ID=CMCHYDRO) 

GETPF (TAPEI,  MIXEDPHASEBOUNDARY,  ID=CMCEOS) 

LIBRARY  (LIBCMC) 

EDITLIB . 


UPDATE (F , P=A ,N=B ) 

FTN (I =CCMPILE , L=0 , PL=55000 , LCM=I ) 
RFL (160000) 

LGO . 

CATALOG (TAPE7 ,CUHEMI , ID=CMCHYDRO) 
EXIT. 

*ECR 

LISTLIB (*, LIBCMC) 

*EOR 

*EOR 


SAMPLE  PROELEM  2:COPPER  HEMISPHERE  LOADED  WITH  COMP  P  UNCONFINED 
$OPTNS  NPRTOP=2 ,  TIMMAX=900,  IFLGST=1 ,  IPMADJ=1,  ITPHSE=1 


$START 

PK(1)=1, 

PK  (2 ) =0 , 

PK  (3) =-0$ 

$RUN 

TSTOP=  6.0E-06 

,  PRDELT=  2.0E- 

-06, 

KUNITW=7 , 

KUNITR=7 , 

NOSLIP=0 , 

LVISC=1 , 

NMAT=2 , 

NMXCLS=800 , 

NTRACR=5 , 

NTPMX=90  0 ,  REZ  =  0 

CYCPH3=-1 , 

NLINER=1 , 

NSLD=300 , 

MAPS=1 , 

IMAX=60 , 

JMAX=200 , 

DMIN=1 00 , 

IPR=100 , 
NFRELP=1 00, 

NDUMP7=1, 

IPCYCL=0 , 

11=3,  12=53, 

CVIS=-1 , 

IGM=0$ 

$GRID 

DX=60* . 05 ,  DY= 

200*. 05$ 

$MTRL 

MFLAG  (2) =0 , 

$ 


MAT ( 1 ) =  3,  RHOIN  (I )=8 .94 , 

MAT ( 2 ) =  121,  RHOIN (2) =1 . 717 , 

CZERC (l)=2.35E+09,  STK1  (1 ) =6 . 95E+10 ,  STK2 (1 ) =5 . 50E+10 , 
STEZ(l)  =5 . 30E+09 ,  RMUI(l)  =4.55E+11,  AMDM(1)=  .9785, 

AMDM (2) =1 .0$ 
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$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GECM 

$GECM 

$GEOM 

$GEOM 

$SLIP 


$DETN 

$DETN 

$DETN 

$DETN 

$DETN 


LT=4  ,  MPN=1 ,  NPTS=300 , 

TH1=3. 14159,  TH2=1. 57080, 

SLA=0 . 0 ,  .  SLB=6 . 35 ,  XR=1.905$ 

LT=1 ,  MPN=1,  NPTS=50 , 

X1S=1 . 905 ,  Y1S=6 . 35 ,  X2F=1.7907, 

LT=-4 ,  MPN=1,  NPTS=300 , 

TH1=1. 57080,  TH2=3. 14159, 

SLA=0 . 0 ,  SLB=6 . 35 ,  XR=1.7907$ 

LT=1 ,  MPN=2 ,  NPTS=100 , 

X1S=0.0,  YlS=2 . 54 ,  X2F=2 . 159 , 

LT=2 ,  MPN=2 ,  NPTS=200 , 

XlS=2 . 159 ,  Y1S=2 . 54 ,  X2F=2.159, 

LT=1 ,  MPN=2 ,  NPTS=50 , 

X1S  =  2 . 159 ,  Y1S=6 . 35 ,  X2F=1.905, 

LT=-4 ,  MPN=2 ,  NPTS=300 , 

TH1=1. 57080,  TH2=3. 14159, 

SLA=0 . 0 ,  SLE=6 . 35 ,  XR=1.905$ 

LT=4 ,  MPN=3 ,  NPTS=300 , 

TH1=3. 14159,  TH2=1. 57080, 

S  L A  =  0 . 0 ,  SLB=6 . 35 ,  XR=1.7907$ 

LT=1 ,  MPN=3,  NPTE=50 , 


Y2F=6 . 35$ 


Y2F=2 . 54$ 
Y2F=6 . 35$ 
Y2F=6. 35$ 


XlS=l . 7907 ,  YlS 
LT=1 ,  MPN= 

XlS=l . 905 ,  YlS 
LT=2 ,  MPN= 

XlS  =  2 . 159  ,  YlS 
LT=-1 ,  MPN= 

X1S=2 . 159 ,  YlS 
LT=1 00$ 

MASTRD  (1)=1 , 
NBGSD(l)  =0, 
MASTRD  (2) =0 , 
NEGSD  (2) =301 , 
IDET=1 ,  NDETS= 
IDET=2 ,  NDETS= 
IDET=0  $ 


=6.35,  X2F=1 

3,  NPTS=50  , 
=6.35,  X2F=2 

3,  NPTS=200 , 
=6.35,  X2F=2 

3,  NPTS=100 , 
=2.54,  X2F=0 

NSLAVD (1)=0, 
NENDMD (1 ) =300 , 
NSLAVD (2) =2 , 
NENDMD (2) =0 , 

2,  XDETS=0 . 0 , 
2,  XDETS=1 . 82 , 


.  905, 
.  159, 
.159, 
.0, 


Y2F=6. 35$ 
Y2F=6 . 35$ 
Y2F=2. 54$ 
Y2F=2 . 54$ 


NBGMD (1 ) =1 , 

NENDSD ( 1 ) =0 , 

NBGMD  (2) =0 , 

NENDSD (2) =600$ 
YDETS=2 . 54 ,  TIME=0.0$ 
YDETS=5 . 62 ,  TIME=0.0$ 


XMIN=0 . 0 ,  YMIN=2 . 54 , 
XMIN=0 . 0 ,  YMIN=5 . 59 , 


XMAX=2 . 159 ,  YMAX=6 .35$ 
XMAX=2 .159,  YMAX=6 .35$ 
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APPENDIX  C:  SAMPLE  OUTPUT 
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10/06/80  SCOPE  2.1.5  8  R  L  VEP  Oil  *♦*  11/06/80  00311 

STS  DEVICES  110/  4/PF  FLS-377K  FU-1750K  MXS-300K  MXL-1305K  MXB-13058 


HH.ttfl.SS  CP’J  SECOND  ORIGIN 
08.2n28.MFA.  ** 
16.08.32  C0000.003  MFZ  • 

16.06.32  00003.006  JOB • 

16.08.33  00030.028  JOB. 

16.00.33  00000.032  JOB. 

16.08.33  00000.033  JIB. 

16.06.33  00003.034  JOB. 

16.08.33  00000.035  JOB. 

16.08.33  00000.03”  MF  Z • 

16.08.33  000  30  .003 
16.08.33  00000.004 
16.08.34  00000.02” 
16.03.34  00030.033 
16.08.34  00000.036 
16.08.34  00030.038 
16. 00.34  00030.052 
16.03. 35  00000.041  LOD. 

16.09.37  00003.059  USR. 

16.09.  37  00000. 0”3  USR. 

16.03.37  00000.073  US*. 

16.09.3”  03000.  O'*  4  USR. 

16.09,37  03330. 0”4  L  HD. 

16.09.37  03000.088  MFZ. 

16.09.36  00030.003 
16.09.36  00000.004 
16.39.4C  30030.027 
16.09.40  30030.033 
16.09.40  00000.036 
16.09.40  00000.038 
16.09,42  00000.087 
16.10.  43  033  34.090  tJSR. 

16.10.43  00004. 0°9  L3D. 

16.17.47  00029.455  USR. 

16.17.47  00029,456  JOB, 

16.1”. 50  03029.456  LID. 

16.25.48  00329.963  MFZ. 

16.25.48  30029.966  MFZ , 

16.25.40  00029.966  USR. 

16.26.00  00040.586  MFZ. 

16.26.00  00000.003  MFZ 
16.26.00  33030.004  JOB , 
16.26.02  00000.02”  103. 
16.26.02  00000.033  UFZ. 
16.26.32  00030.036  UFZ. 


BP.L  HUS/ BE  1.3  L 499  VEP.  007 
—LCTI L#  STUFZ#  T1000#  PO. 

-ACCTUUT(MD***I 
“REQUEST! TAPE?#  *PF) 

-GET PF !  LIBCMC#  LIBCMC#  ID-CHCEOS) 

-GET°F  !  A#  BP.LHEL 0#  ID-C  MCMYDRO) 

-GETPF ( T  APE1#  MIXED PHASE  BOUNDARY#  ID-CMCEOS) 

— L IBRARYCL IB CMC ) 

JM262  -  STAGE  PF  ST-MFZ  LFN-LIBCHC 


10/00/80 

SAMPLE  PROBLEM  It  CU  HEMISPHERE  LACETERA 


MFZ 
JOB. 

LOO . 

MFZ. 

UFZ . 

MFZ. 

MFZ. 

-COITLIB. 

FORTRAN  LIBRARY  508 
STOP 

042100  FINAL  EXECUTION  FI. 

.014  CP  SECONDS  EXECUTION  TIME. 
-UPDATECF  /  ”*A#fUB) 

JM262  -  STAGE  PF  ST-MFZ  LFN-A 


-SJACP#STMFZ#T100#P?000.  ATT/CAT/PRG  SPOT  JOB 
—ACCOUNT  (CS***)  SYSTEM  JOB  5JACP 
-PFCONV.  CALI  ATTACM/CAT/PURGE  UTILITY 

PF646  -  PTMACRD  -  ATTACH  -  LIBCMC  -  LIBCMC 
PF254  -  CYCLE  10  ATTACHED  FROM  SN-SYSTEM 
JM262  -  STAGE  PF  SPOT  LFN-PF76C 

JM262  -  STAGE  PF  SPOT  LFH-PF76C 


05/27/80 


MFZ. 
JOB. 
L  Tj. 
MFZ. 
MFZ. 
MFZ. 
UFZ. 


-SJACP#STNFZ,T100#P7000.  ATT/CAT/PRG  SPOT  JOB 
-ACC  OUMT ( C  S*  ** )  SYSTEM  JOB  SJACP 
-PFCONV.  CALL  ATTACH /C AT/PUP GE  UTILITY 

PF646  -  PFMACRf)  -  ATTACH  -  A  -  BRLHELP 

PF254  -  CYCLF  23  ATTACHED  FROM  SM-SYSTEM 
JM262  -  STAGE  PF  SPOT  LFN-Pf 76C 

JN262  -  STAGr  PF  SPOT  LFN-PF76C 

UPDATE  COMPLETED 

“FTN! I -COMPILE #L»0#PL* 55000# LCM- I ) 

25.34”  CP  SECONDS  COMPILATION  TIME 
-P.FL  ( 160300) 

-ICO. 

LD610  -  FLS  REQUIRED  TO  LOAD  -  0043770  OU.COG 
LD603  -  EXECUTION  INITIATED  OS. EXP 
FORTRAN  LIBRARY  510  09/30/80 

J  1262  -  STAGE  PF  ST-MFZ  LFN-TAPE1 

-SJACP#  ST  MFZ#  710  0#  P700Q.  ATT/CAT/P  P.G  SPOT  JOB 

-ACC  UUMT ( C  5*  ** )  SYSTEM  JOB  SJACP 
-PFCONV,  CALL  ATTACH /C  AT/ ftUP  GE  UTILITY 

PF646  -  PF MACRO  -  ATTACH  -  TA*E1  -  MIXEDPHAS EBQUMO ART 
PF254  -  CYCLE  19  ATTACHED  FROM  SN-SYSTEM 


16,26.02 

03000 

.036  MFZ. 

JM262  -  STAGE  PF  SPOT 

LFN-PF76C 

16.26.02 

00030 

.045  MFZ . 

JM262  -  STAGE  PF  SPOT 

LFN-PF76C 

17.34.06 

03897.272 

’JSP. 

EXIT 

17.34,06 

00897,272 

USR . 

1690  ^0  FI  UAL  EXECUTION  FL. 

17.34.06 

008 97. 2**3 

USR. 

967.305  CP  SECONDS  EXECUTION  TINE. 

17.34. 06 

00397.276 

JOB. 

-CATALOG!  TAP  E7#CUUEM  I,  ID-C  HCMYDP.O) 

17.34.06 

00897.281 

MFZ. 

PF060 

-  CYCLE  3  CATALOGED  ON 

SN-SYSTEM 

17.34.06 

03097.282 

JOB. 

-EXIT. 

17.34.06 

00897.204 

UFZ. 

JU166 

-  MAXIMUM  USER  SC M 

1600008 

WORDS 

17.34.06 

09C9n285 

MFZ. 

JM16* 

-  .MAXIMUM  USER  LCM 

651000B 

WOROS 

17.34.06 

00897.205 

MFZ. 

J  Ml  70 

-  MAXIMUM  JS+IO  LCM 

2238 

BUFFERS 

17.34.06 

00097,205 

MFZ. 

P.U””0 

-  MAXIMUM  ACTIVE  FILES 

8 

17.34.06 

00897.265 

UFZ. 

F,  1771 

-  OPEN/CLUTE  CALLS 

48 

17.34.06 

00397, 285 

MFZ. 

RU772 

-  DATA  TRANSFER  CALLS 

193#  876 

17,34.06 

00097. 2B5 

MFZ. 

RM7”3 

-  CONTP OL /PDS ITIDNING  CALLS 

1#  002 

tT.it.n/. 

mt  7 . 

-  ti h  r\KTk  TDAu^rrp  mi* 

1 n. 066 
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sample  problem  2  -  copper  hemisphere  loaoeo  with  comp  b  me  umconfimed 
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detonation  time  CALCULATION  FOP  EXPLOSIVE  PACKAGES 
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PACKAGE  *m'lAl  INITIAL  CONDITIONS 

NUMBER  DENSITY  DENSITY  S.t.E.  V  V  MATERIAL 

(RHOZI  (RM3IN) 
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PRESSURE  CYCLE*  69.0  TIME -2 .OOOOOE-Q6  SECONDS 
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